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fragments were too low to induce a ground-shock condition suffi- 
cient for cratering. However* we do suggest that this body developed 
bow-shock peak overpressures and shock-initiated winds in the 
atmosphere near the ground sufficient to cause incipient distur- 
bances on the venusian surface, depending on rheologic states of the 
target materials. 

The 10 km-diameter body traveled from 45 km altitude down to 
the surface in 2 s with less than 2% deceleration and would have 
caused major atmospheric, impact-cratering, and ejecta effects. We 
estimate that the final crater would have been 1 50 to 200 km across 
but only a few kilometers deep, which we infer would be due in part 
to the high temperature, low strength, rapid relaxation rate of the 
target rock. This large body generated an enormous bow-shock 
wave with pressures in excess of 600 kbar that reached the ground 
surface and swept outward at over 10 km/s. It also evacuated a 
region over 40 km in diameter along the trajectory of the body in the 
atmosphere. This hot (-20,000-200,000 K), low-density region 
continued to expand and remained open for over several minutes 
after impact with substantial amounts of high-angle ejecta traveling 
through it to high altitudes. The effects of the 10-km body in Venus’ 
environment are still under analysis, but several of the major aspects 
appear similar to the asteroid-impact results calculated for the 
Earth. The energetics of such large events tend to initially over- 
whelm both atmospheres. 

In summary, our numerical simulations indicate that for weak, 
cometlike objects impacting into Venus’ atmosphere at 20 km 7s, a 
100-m -diameter body completely disintegrates at high altitude and 
produces no ground surface effects. A 1000-m -diameter body also 
disintegrates but produces a strong bow-shock wave immediately 
followed by a dispersed band of vapor, melt (?), and fragments that 
reach the ground to produce strong shock pressures and radial winds 
in the near-ground atmosphere. We suggest that these conditions 
could induce incipient surface disruption and produce subtle surface 
features, depending on the rheologic states of the target materials. 
We suspect that the dark and light splotches that are several tens of 
kilometers in diameter seen in the Magellan SAR images [3-8], and 
believed to be impact induced, could have been generated by 
cometlike bodies a few kilometers in diameter. The bodies could be 
somewhat smaller if the splotches were produced by stronger 
asteroidlike bodies that fragmented lower in the atmosphere; these 
cases remain to be explored. 

References: [ 1 ) Seiff A. (1 98 3) In Venus, 1 045-1048, Uni v . of 

Arizona, Tucson. [2] Roddy D. J. et al. (1987)M/.7. Impact Engng., 
525-541. [3] Arvidson R. E. et al. (1991) Science , 270-275. 
[4] Phillips R. J. et al. (1991) Science, 288-297. [5] Zahnle K. J. 
(1991) Eos, 44, 289. [6] Zahnle K. J. (1992) LPSC XXIII, 1565. 
[7] Soderblom L. A. et al. (1992) LPSC XXIII, 1329-1330. [8] 
Schaber G. G. et al., in press. \ 

N93-14366 Ms M/i 

THE EFFECTS OF VENUS’ THERMAL STRUCTURE ON 
BUOYANT MAGMA ASCENT. S. E. H. Sakimoto and M. T. 
Zuber, Department of Earth and Planetary Sciences, The Johns 
Hopkins University, Baltimore MD 21228, USA. 

The recent Magellan images have revealed a broad spatial 
distribution of surface volcanism on Venus [1]. Previous work in 
modeling the ascent of magma on both Venus and Earth [2-5] has 
indicated that the planetary thermal structure significantly influ- 
ences the magmatic cooling rates and thus the amount of magma that 
can be transported to the surface before solidification. In order lo 
understand which aspects of the thermal structure have the greatest 


influence on the cooling of ascending magma, we have constructed 
magma cooling curves for both plutonic and crack buoyant ascent 
mechanisms, and evaluated the curves for variations in the planetary 
mantle temperature, thermal gradient curvature with depth, surface 
temperature gradient, and surface temperature. The planetary ther- 
mal structure is modeled as 



( 1 ) 


where T is the temperature, T 0 is the source depth temperature, x = 
1- (Tj/T 0 ) where T s is the planetary surface temperature, Z is the 
depth, Zo is the source depth, and n is a constant that controls thermal 
gradient curvature with depth. Equation (1) is used both for math- 
ematical convenience and flexibility, as well as its fit to the thermal 
gradients predicted by the cooling half-space models [6]. Wc 
assume a constant velocity buoyant ascent, body-averaged magma 
temperatures and properties, an initially crystal-free magma, and 
the same liquidus and solidus for both Venus and Earth. 

The cooling model for the plutonic ascent has been described in 
detail in earlier publications [2-5], and is a low Reynolds number, 
high Peclet number problem of heat transfer through a thin thermal 
boundary layer around a sphere. The resulting plutonic cooling 
curves, which are dominated by the convective cooling terms and 
strongly influenced by the planetary thermal structure, are then 
expressed mathematically by 
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Fig. 1. Example of a typical cooling curve plot forpluton ascent. The curves 
are plotted for a given thermal structure and/or planet, and are contoured with 
their associated dimensionless ascent (J^) values. In order for the magma to 
reach the surface unsolidified, the cooling curve must not cross the solidus. 
This plot is for a dry olivine tholeiitc. The liquidus and solidus are obtained 
from [11}. Thermal structure parameters: n = 2; dT/dZ= i .0 (dT/dZof Earth); 
Zq= \ .0 (Zo of Earth). 
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where J = 3NuK/a 2 , y = agV/C pi t is time, 1$ is total ascent time, 
Nu = 0.8Pe 1/2 , Pe = Va/K, Pe is the Peclet number, V is the velocity 
of magmatic ascent, a is the body radius, K(=l x 1(H mV 1 ) is the 
thermal diffusivity , a (=6 x IQ -5 deg -1 ) is the coefficient of thermal 
expansion, g is the gravitational acceleration, and Cp(=1.25 x 1CH 
ergs Kg-^K" 1 ) is specific heat capacity. T is the mean magma 
temperature, T 0 is the magma temperature in the source region, and 
n is a constant that defines the shape of the planetary thermal 
gradient (equation (1)). Equation (2) reduces to the expression for 
the thermal gradient (equation (1» for an infinitely slow ascent 
(dimensionless ascent time Jto = ®°), and to the adiabatic curve 
T/T 0 = e - ? 1 for an infinitely fast ascent (Jto = 0). A typical plot of the 
resulting cooling curves for terrestrial conditions, contoured in Jto 
values, is illustrated in Fig. 1 . In order for the magma to reach the 
surface un solidified, the cooling curve must not cross the solidus 
before it reaches the surface. The allowable Jto values obtained from 
the cooling curve plots for Venus and Earth can be directly com- 
pared to obtain relative minimum magma ascent velocities, source 
depths, and body sizes. The results are shown in Fig. 2. 

The cooling model for the buoyant crack ascent has previously 
been described briefly in [2]. It is the problem of a magma at an 
initial temperature T 0 placed in contact with the wall rock of 
temperature T m . This problem was initially solved by [7], and their 
solution for the average temperature T as a function of time is 
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where the notation is the same as in equation (2), and the right side 
is constant fcrr any single dimensionless ascent time (Kt/a 2 ). This 
result is for a constant wall rock temperature, but can be adapted to 
a variable wall rock (thermal gradient) temperature by approximat- 
ing an incremental magma ascent in a simple numerical scheme 
where the initial magma temperature T 0 at any location m is, instead 
of the source depth temperature, the final magma temperature at the 
previous location m-1 [2]. If heat is conducted ahead of the magma 
body , the boundary temperature of the magma and wall rock will not 
be constant (T^^t = T m ) at any given location, but will be the 
average of the two initial temperatures of the magma and wall rock 
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Fig. 2. Plot of the results for basaltic pluton ascent on Venus. The results are 
nondimensionalized by the corresponding terrestrial values, and the source 
depth values are the direct result of the assumed planetary thermal structure 
(equation (1)) . For example, for a Venus thermal gradient ihe same as the 
terrestrial thermal gradient and n = 2, the corresponding minimum source 
depth is about 60% of the terrestrial value the minimum ascent velocity is 
about 30% of Earth’s and the minimum body size is about 70% of Earth’s. 
Thermal structure parameters: n = 2; dT/dZ = 1.0 (dT/dZ of Earth); To = 1*0 
(Zo of Earth). 


at any location for the majority of the cooling time [8]. For this case, 
the contact temperature at the mth position is the average of the final 
temperature of the magma at the m— 1 position (Tjj^i) and the initial 
wall rock temperature at the mth position (T^). The preliminary 
results from this model indicate that the effect of the planetary 
thermal structure is of the same order of the effect seen in the pluton 
model. 

In general, for both ascent mechanism models presented here, 
the influence of the planetary thermal structure parameters for 
Venus in the probable order of decreasing importance is surface 
temperature, surface temperature gradient, thermal gradient curva- 
ture with depth, and planetary mantle temperature. The higher 
surface temperature of Venus, for otherwise similar planetary 
thermal structures, allows considerably smaller minimum possible 
crack sizes and/or magma body sizes, and slower ascent velocities 
than would be possible on Earth for a reasonable range of Venus 
source depths and surface thermal gradients. This surface tempera- 
ture effect is greater for more primitive magma compositions, and 
may be greater for magmas of higher crystallinity. A higher venu- 
sian surface thermal gradient has the same effect of the higher 
surface temperature on magma transport, but to a much lesser 
degree. Similarly, for higher values of thermal gradient curvature 
with depth (higher n in equation (1)), the minimum possible ascent 
velocity and body/crack size also decreases slightly. If the mantle 
temperature for Venus is elevated by a hundred degrees or so over 
that of Earth [9], it should result in a modest increase of melt 
production and magma transport to the surface compared to Earth. 
The effectof the range of Venus surface temperatures with elevation 
(390°-470°C or 660°-740°K [10]) is under investigation, and is also 
anticipated to have a significant effect on magma transport, possibly 
greater than that of the higher mantle temperature. 
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Though there is no plate tectonics per se on Venus [1], recent 
Magellan radar images [2] and topographic profiles [3] of the planet 
suggest the occurrence of the plate tectonic processes of lithospheric 
subduction [4] and back-arc spreading [5] . The perimeters of several 
large coronae (e.g., Latona, Artemis, and Eitliinoha) resemble Earth 
subduction zones in both their planform and topographic profile. 
McKenzie et al. [4] have compared the planform of arcuate struc- 
tures in Eastern Aphrodite with subduction zones of the East Indies. 
The venusian structures have radii of curvature that are similar to 
those of terrestrial subduction zones. Moreover, the topography of 
the venusian ridgc/trench structures is highly asymmetric with a 
ridge On the concave side and a trough on the convex side; Earth 
subduction zones generally display this same asymmetry. 




